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Abstract 

We discuss the usage of continuous external momenta for computing renor- 
malization factors as needed to renormalize operator matrix elements. These 
kind of external momenta are encoded in special boundary conditions for the 
fermion fields. The method allows to compute certain renormalization fac- 
tors on the lattice that would have been very difficult, if not impossible, to 
compute with standard methods. As a result we give the renormalization 
group invariant step scaling function for a twist-2 operator corresponding to 
the average momentum of non-singlet quark densities. 



1 



1 Introduction 



One of the most important contributions lattice gauge theory can provide to test 
QCD and to interpret experimental data, is the computation of non-perturbatively 
renormalized matrix elements. The example we are interested in here are matrix 
elements that are connected to moments of parton distribution functions as can be 
extracted from global fits to experimental data in deep inelastic scattering. The 
matrix elements are related themselves to certain operator expectation values which 
are accessible from lattice simulations. 

When considering operators that contain derivatives, a saturation with external 
momenta is needed to perform the necessary contractions. This applies for comput- 
ing the renormalization constants as well as matrix elements of such operators. On 
an euclidean lattice, the standard momenta are quantized in units of 2ira/L with a 
the lattice spacing and L the linear extent of a lattice of physical size L 4 . Unfor- 
tunately, the introduction of a momentum in a numerical simulation leads often to 
either large lattice artefacts or to large statistical uncertainties or to both such that 
reliable measurements of such quantities become difficult. 

In this paper we will demonstrate how the introduction of special boundary con- 
ditions of the fermion fields, first advocated in P, can help in such a situation. The 
boundary conditions lead to continuous values of external momenta ad/L where 9 
varies continuously in the range < 9 < 2n giving a large flexibility of a momen- 
tum definition in lattice simulations. In particular, the minimal value of a lattice 
momentum can be chosen smaller than in the standard set-up which might lead to 
smaller lattice artefacts. As we will see in the example discussed below, the use of 
these kind of continuous external momenta will allow us to compute renormalization 
constants reliably on the lattice for which measurements were very difficult in the 
standard set-up. It might be, however, that the method is much more general and 
goes beyond the application investigated here. 

The physical example we will discuss in this paper is the calculation of renormal- 
ization constants as needed for twist-2, non-singlet quark operators in the framework 
of computing moments of parton distribution functions in deep inelastic scattering. 
The interest in this calculation is twofold. First, we want to demonstrate the feasi- 
bility of using continuous momenta in a lattice computation. 

Second, the success of this demonstration will allow us to correct a small mis- 
match of our earlier work: In a series of papers (21 El IH El E] we have demonstrated 
that the Schrodinger functional (SF) formalism [7j can be used to compute non- 
perturbatively moments of parton distribution functions (for summaries of the re- 
sults, see [HI EE EI])- A problem that arises naturally in lattice computations of such 
quantities is that for a given continuum operator there is not a unique representation 
of that operator on the lattice. In fact, in the old work cited above, the operator 
used to compute the matrix element and the operator used to compute the corre- 
sponding renormalization constant were in two different representations. Although 
in ref. jH] we argued that this mismatch only leads to a small, negligible error, the 
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situation remains rather unsatisfactory, as we want a fully non-perturbative way to 
evaluate the physical matrix element without systematic uncertainties. The usage 
of the continuous external momenta will allow us to eliminate this small systematic 
uncertainty. In particular, we will show that we can compute the scale dependent 
renormalization constant and the corresponding (ultra-violet) invariant step scaling 
function for a particular lattice representation which is needed to perform the cor- 
rect renormalization for the physical matrix element. In this way, we will be able to 
provide a non-perturbative lattice computation of the average momentum carried 
by the quarks in a hadron. 



2 Basic definitions and choice of boundary condi- 
tions 

A bare, local operator barc as considered here, is renormalized multiplicatively by 
a scale fi dependent renormalization constant Zo(fi). The renormalized operator 
R (fi) is then given by 

O R Qi) = Z?0*)O* m . (1) 

The evolution of O from a renormalization scale fii to (12 is described by the con- 
tinuum step scaling function oz a that can formally be written as 

Zo(// 2 ) 

VZoWW = , r ■ (2) 

The individual Z's in eq. (J2J) are only well-defined within a given regularization 
scheme and would be divergent when the regularization is removed. The step scaling 
function, however, is a well-defined quantity even in this limit. In the following we 
will give a rigorous definition of the renormalization constants and of the step scaling 
function using the lattice regularization in the Schrodinger functional (SF) scheme. 

The SF scheme is based on the formulation of QCD in a finite space-time volume 
of size L 3 ■ T. In this paper we will always use L = T. In this scheme a change of 
renormalization scale amounts to a change of the box size L at fixed bare parameters. 
By considering a sequence of pairs of volumes with sizes L and sL, one can study the 
evolution of a given local operator under repeated changes of the scale by a factor of 
s. Effectively in this way one builds up a non-perturbative renormalization group. 
To be specific, we will from now on work on a euclidean lattice as regulator. 

The fermion and gauge fields on the lattice are defined in the standard way, ful- 
filling SF boundary conditions, as detailed in ref . [7] . Besides these special boundary 
conditions in time direction, we will impose generalized boundary conditions in the 
spatial directions (denoted by k) for the fermion fields ip(x) [I] 

ip(x + Lk) = e i9k ?p(x), ip(x + Lk) = e"^(x), k = 1, 2, 3, (3) 
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while the gauge fields U(x, k) are chosen to be periodic in the space directions. The 
values of 9k can be chosen in the interval < 9k < 2tt. For 9k = we obtain periodic 
and for 9k = 7r we get anti-periodic boundary conditions. 

The generalized boundary conditions of eq. (J3J) can be implemented in the defi- 
nition of the gauge covariant lattice derivatives: 

V '^{x) = ^ [\ tl U{x,ij)ip{x + ap) -ip(x)] 

V*?/>(:r) = - [i/j(x) — \*JJ{x — ap, ^)^(x — ap)\ (4) 



with fi = (k, 4), 



A M = e^/ L , 0<9 k <2n, 9 4 = 0, (5) 



and a the lattice spacing. From a technical point of view the fields ip(x) and ip(x) 
are then implemented with the usual periodic boundary conditions in the space 
directions utilizing this generalized definition of the covariant derivative of eq. (j3J). 

The crucial observation pQ is that the factor e taek / L can be interpreted as an 
external momentum with the intriguing property that it can assume continuous 
values, in contrast to the standard, quantized lattice momenta that assumes values 
in units of 27ra/L only. 

In order to explore the flexibility of a momentum definition given by the gen- 
eralized boundary conditions in eq. (J3J) we have concentrated in this work on the 
twist-2, non-singlet (quark) operator. This amounts to consider operators of the 
form (the flavor structure is specified by the Pauli matrix r 3 ) 

Of, u {x) = -$(x)^ {lJi D u } )^~ z 4>(x) - 5^ ■ trace terms . (6) 

where { • • • } means symmetrization on the Lorentz indices and 

5„=5 M -5 M ; ^ = i[v M + v;]. (7) 

There are two representations of such a non-singlet operator on the lattice jllj . 
The first representation takes /i ^ v whereas the second uses fi = v. The precise 
definitions of the operators used here are 

12 (x) = ^(x)7{i B 2} \r^(x) (8) 



and 

3 



3 , , 



1 r <-» 1 „ l ^ 

O u (x) = -$(x) 74 Da -- 7* D k —ij>(x) . (9) 



In both cases an external momentum has to be supplied to compute their renor- 
malization constants. We will realize these momenta with non-vanishing values of 
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6k, implicitly given through the covariant derivatives in the fermion action and the 
operators in eqs. (jHJ), ©• The precise choice we adopt here is 

9 = Bx + 0, 9 2 = 9 3 = . (10) 

In order to make the discussion in the following self-consistent, let us recall the 
definition of bare correlation functions in the SF scheme [3j as needed to compute 
appropriate renormalization factors. The correlation function of a given operator O 
is given by 

f o {x A /L,0) = — ^(o^x) C(y)r^r 3 c(z)) , (11) 

x,y,z 

where T is a Dirac matrix that in our case is 72 for the O12 operator, while it is 71 
for (9 44 . In eq. (fTTj) ( and ( are classical boundary fields at 24 = and we sum over 
space. In order to normalize such correlation functions properly, we also need the 
"boundary to boundary" correlation function fx, 

= £ <C>)75^rr(v)C(y)75^C(z)> . (12) 

u,v,y,z 

Note that the classical boundary source fields (, ( at X4 = and (' at X4 = L are 
renormalized multiplicatively with a common renormalization constant [T2*| HHj. 
The momentum (9) dependence of the correlation functions fo in eq. (fTTJ) and f\ in 
eq. ()12j) is again implicitly given through the definition of the covariant derivative 
eq. that appears in the fermion action and in the operator considered. 

It is instructive to look at the 6*-dependence of the correlation functions fo 44 and 
fo 12 for fixed values of X4/L = 1/2 and X4/L = 1/4 already at tree- level. Fig. [T] 
shows the tree-level correlation functions fffl. It can be observed that the signal 
strongly depends on the chosen value of 9, the definition of the operator and the 
choice of x&/L. For small values of 9 all correlation functions show a linear behavior 
in 9 and vanish at 9 = 0. In contrast to the correlation function of the operator 
O12, the correlator for C 44 , which does not depend on x^/L at tree-level, becomes 
very small for values of, say, 9 > 4.5, see the inlay in fig. ^ This implies that for 
a real simulation there is a danger that for these ^-values the signal may become 
very small, too. In addition, it seems that for these values of 9, due to large lattice 
artifacts (cf. sect. 2.2) it may become hard to extract even the one-loop anomalous 
dimension. It should be concluded therefore that the choice of 9 in a simulation needs 
an optimization procedure. Before we can start a discussion of how to perform this 
optimization we first want to give here the renormalization conditions we have used 
in order to determine the renormalization factors which are defined as follows. 

Z {a/L,xjL,9,9') = c f (a/L, x 4 /L, 9)/ y/f^a/L, 9') 
Z {a/L,x 4 /L,9) = c f {a/L,X4/L,9) 
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Figure 1: ^-dependence of fo at tree- level. For the O44 representation the tree- level 
values at x 4 = L/2 and x 4 = L/A are the same. 

where the correlation functions are being evaluated at vanishing quark mass m q = 
and for vanishing boundary gauge and fermion fields. The renormalization condi- 
tions are chosen such that Zq = 1 and Zq = 1 at tree-level of perturbation theory 
at a fixed scale fi = 1/L [21 IS] while the external parameters 9, 9' and X4/L can be 
varied. From these conditions we obtain 



f (o) 



AO) ' " f (o) • 
Jo Jo 

The quantities c and c are the expressions of the correlation functions in eq. (jllj) 
and in eq. (|12|) at tree-level with the corresponding arguments. Note that the actual 
renormalization constant needed is Zq whereas Zq is only an auxiliary quantity 
that, however, proved useful in our previous work 0. Note that in the definition 
of the renormalization constants in eq. (fTHj) we leave the freedom to choose 9^9' 
for fo and fi. It is important to stress here that a different choice of x$/L, 9 and 
9' defines a different renormalization scheme. In this paper we have considered two 
families of renormalization schemes, 

• scheme A 9' = 9 

• scheme B 9' = . 
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Both schemes are studied in the following with respect to analyzing the convergence 
of perturbation theory, the size of cutoff effects and the signal to noise ratio. In these 
studies the external parameters 9 and x±/L are tuned to optimize these criteria. 

2.1 Convergence of perturbation theory 

In this section we present a perturbative analysis of renormalization constants for 
the local operators in eq. (JEJ) and in eq. In particular we compute the 2-loop 
anomalous dimension for the Z-factors for the two schemes discussed above. In bare 
perturbation theory the renormalization constants have an expansion of the form 

oo 

Z o (g , a/L, Xi /L, 9) = 1 + J2 Z ( k) (a/L, xjL, 9) g 2 k , (15) 

k=l 

where in the limit a/L — > the coefficients Z$ are polynomials in ln(L/a) of 
degree k up to corrections of 0(a/L). In particular the coefficient of the logarithmic 
divergence in Zq is given by the one-loop anomalous dimension j , and thus Z$ 
is parameterized as 

Zg ] = B (8, x A /L) - 7o ln(L/a) + 0(a/L) . (16) 

The values for the 1-loop anomalous dimension 70 and the constant piece Bq{9, X4/L) 
can be computed in lattice perturbation theory for the SF scheme following the tech- 
niques developed in [21 E]- For the matching of the non-perturbative lattice results 
and perturbation theory at very high energies, it is important to know also the 2- 
loop anomalous dimension ji*. Since 71 is not universal, it becomes necessary to 
compute the dependence of 71 on the values of 9 and x±/L. In this way it becomes 
possible to control for which choice of these parameters 71/70 is small in order to 
have a good behavior of the perturbative series. 

The coefficient Z a oc 7iln(L/a) can actually be obtained without an explicit 
2-loop calculation in the SF scheme if the 2-loop anomalous dimension is already 
known from a different renormalization scheme ^H] • The formula relating the 2-loop 
anomalous dimensions in the SF and the MS schemes is given by 

7l = 7 fS + 2bo AZ 0) _ 1qX (1) . (17) 

In eq. ()17j) 60 is the universal 1-loop coefficient of the /3-function, AZq^ is the 1-loop 
difference of the renormalization constants from a finite renormalization that relates 
two mass-independent schemes and X% is the perturbative factor that relates the 
renormalized couplings in the two schemes considered (see the appendix for explicit 
expressions of these quantities). 

A subtlety is that the factor AZ^ needs to be computed in two different regu- 
larizations, in order to avoid a calculation with the SF scheme using a dimensional 
regularization (DR). 

*When we omit an explicit scheme index, we always mean the SF renormalization scheme. 
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Renormalization 
Lattice reg. schemes Dimensional reg. 



SF 



MS 



MOM 



t 



Figure 2: Graphical representation of the finite renormalization between schemes 
when a different regularization is used. 

Let us explain this fact. A matrix element O s R renormalized in a certain scheme 
s is obtained by 



O 



R 



rrs,reg srpeg 



where O r b e9 is the bare matrix element computed within a certain regularization 
reg and Z s ^ e9 is the renormalization constant that depends on the renormalization 
scheme s used and on the regularization reg. Operators renormalized in two different 
schemes but using the same regularization can be related by a finite renormalization 



O s AZ, 



o 



(19) 



One important observation here is that AZ S ^ S is independent from the regular- 
ization used to compute the renormalized matrix element and the corresponding 
renormalization constant. This allows to relate two schemes: in principle it is possi- 
ble to compute the complete 1-loop renormalization constant (anomalous dimension 
and finite part) in the SF scheme, using the dimensional regularization and then 
connect directly with the MS scheme (dashed line in fig. |2J), 



AZ§T 



-MS 



J o 



l /Z { 



(20) 



A lattice perturbative computation will give, on the other side, important informa- 
tion about the amount of discretization errors at 1-loop and would help the numerical 



7 



9 


A / n\ 

7i (0) 


7i (0) 


7i (0)/7o 


7i (0)/7o 


0.1 


06584(1") 


06551(1) 


1 4621(6) 


1 4548 ( 3) 


0.4 


06083(1) 


05669(1) 


1 3509(6) 


1 2590(2) 


0.7 


0.05196(1) 


0.04135(1) 


1.1539(6) 


0.9182(2) 


1.0 


0.04211(1) 


0.02482(1) 


0.9352(6) 


0.5512(2) 


1.3 


0.03313(1) 


0.01029(1) 


0.7357(6) 


0.2286(2) 


1.6 


0.02561(1) 


-0.00139(1) 


0.5687(6) 


-0.0308(3) 


1.9 


0.01951(1) 


-0.01051(1) 


0.4332(6) 


-0.2334(3) 


2.2 


0.01449(1) 


-0.01775(1) 


0.3217(6) 


-0.3943(3) 


2.5 


0.01019(1) 


-0.02375(1) 


0.2263(6) 


-0.5273(3) 


2.8 


0.00627(3) 


-0.02903(3) 


0.1392(6) 


-0.6446(6) 


3.1 


0.00231(3) 


-0.03412(3) 


0.0512(6) 


-0.7578(6) 



Table 1: 2- loop anomalous dimension for the O44 representation computed in 2 different renor- 
malization schemes (cf. text) for several values of 9 and X4/L = 1/2. For X4/L = 1/4 the pattern 
is very similar. 



calculation. In order to avoid an additional computation in the SF scheme one may 
now use the MOM scheme, where a complete 1-loop computation has been done for 
the renormalization constant in both regularization schemes |16j . 

The desired factor AZ§, F *~ m relating the SF to the MS scheme (full line in fig. EJ) 
can now be obtained from 

AZ 5F^MS = AZ SF^MOM AZ MOM^MS ^ 

The two factors AZ are computed using different regularizations, namely AZq F ^ moa 
is computed on the lattice and AZq OM *~ ms in dimensional regularization. Using 
the complete one loop result ^E] i n the MOM scheme which exists in both, lattice 
and dimensional, regularizations, it is then possible to compute Azffi and from this 
finally 71. 

In tabled we give values for 71 for selected values of 9 and x^/L = 1/2. From 
this table we observe that the ratio 71/70 becomes smaller for increasing values of 
9 if we choose the scheme A {6 = 9'). For scheme B (9' = 0), we find a kind of 
minimum of this ratio for 9 « 1.6. We repeated the analysis for the anomalous 
dimensions for X4/L =1/4 and found a very similar pattern. Thus we have a first 
indication for the choice of 9 to be chosen. 



2.2 One-loop cut-off effects 

As discussed above, an essential element to obtain the non-perturbative scale de- 
pendence of a matrix element is the step scaling function of eq. (0) which describes 
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a change from a scale L 1 to (sL) 1 , with a scale factor s. Employing a lattice 
regularization, we define the lattice step scaling function t\ 



-2 it u Z o (g ,sL/a;x i /L, 



E Zo (s,5f ,a/L;x A /L,9) 



Zo(9o,L/a;x 4 /L, 



(22) 

t/(L)=fixed 



where #(£) is the renormalized coupling at scale L 1 . The desired continuum step 
scaling is obtained from a limit procedure 

lim S Zo (a, # 2 , a/ L; x 4 / L, 0) = a Zo (s, # 2 ; ar 4 / L, 0) . (23) 

In the following, the above step scaling functions and the renormalized coupling are 
computed in the SF scheme and the scale factor is set to s — 2. To one-loop order 
of perturbation theory we find 

S Zo (2, f, a/L; x 4 /L, 9) = 1 + k(a/L; x±/L, 9)f + 0(f) (24) 

with 

k(a/L; x 4 /L, 9) = z£ ] (a/2L; x 4 /L, 9) - Z% ] (a/L; x 4 /L, 9) . (25) 



In order to see how fast the continuum limit in eq. (j23j) is approached, we define the 
normalized deviation from this value: 

r / it i t a\ k(a/L;x 4 /L,9) - fc(0) 

6 k (a/L; x A /L, 9) = -— . (26) 

Here k(0; x 4 j L,9) = — 70 hi(2) is the continuum limit value which is independent 
from X4/L and 9. The quantity <5k in eq. (|2l)|) contains all the lattice artifacts at 
0(g 2 ). The results for 44 with x 4 /L = 1/2 are displayed in fig. El Note that 
we have worked in the unimproved theory where we expect the lattice artifacts to 
decrease asymptotically with a rate proportional to a/L. We show our results for 
5k(a/L; X4/L, 9) for various values of 9. The full symbols denote results obtained in 
the renormalization scheme A, while the open symbols correspond to results in the 
scheme B. 

We see from fig. El that there is no real preference for choosing scheme A or 
scheme B as far as the lattice artefacts are concerned. On the other hand, for 
large values of 9, the lattice artefacts become very strong which would lead to 
a difficult continuum extrapolation in the numerical simulations. It seems from 
this perturbative analysis that values around 9=1 should be a good choice for a 
fast convergence of perturbation theory and for keeping lattice artefacts well under 
control. 



2.3 Signal to noise ratio 

The analysis of the tree-level correlation function of the operator (9 44 demonstrates 
that for 9 = 2it, corresponding to the smallest value of standard quantized momenta 
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Figure 3: Discretization errors in the unimproved step scaling function at 1-loop 
order of perturbation theory. Points are for L/a = 6, 8, . . . , 24 in steps of two. Full 
symbols denote scheme A, while the open symbols denote scheme B. 
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Figure 4: ^-dependence of AZ/Z. The chosen value of k corresponds to m q 
The lines only connect the points to guide the eye. 



on the lattice, the value of the correlation function becomes tiny. If this also applies 
in the interacting case, it may become very difficult to obtain a reliable measure- 
ment from a numerical simulation. On the other hand, fig. Q shows that employing 
the boundary conditions of eq. J3J) and choosing 9 < 2n the situation improves 
considerably. 

Of course, the tree-level analysis can not give any information on the ^-dependence 
of the relative statistical error AZ/Z of a renormalization constant, which is the most 
important benchmark for a practical numerical simulation. In order to compute 
AZ/Z for Zq A4 and Zq 12 , we performed therefore a simulation for non-perturbatively 
0(a)-improved (clover) fermions ^7] at /3 = 7.0203, k = 0.134707 (corresponding to 
k c ) on a 16 4 lattice. We show in fig. E]the result of this investigation. 

Let us start the discussion on the relative error for Z 0l2 . For 9 > 1, the relative 
error appears to be almost constant up to 9 = 2tt and independent of choosing 
x^/L = 1/2 or X4/L = 1/4. This means that for the operator 12 choosing a 
physical momentum (as done in our earlier work) or introducing 9 would have led 
to the same quality of the final results. 

For the operator O44 the situation is different, however. At 9 = In the relative 
error of Zq 44 evaluated at x 4 /L = 1/2 is a factor 5-6 larger than the corresponding 
relative error of Zq 12 . For x 4 /L = 1/4 the situation is better, but the error is 
again about a factor 2 worse than for Zq 12 at least for 9 = 2tc. We observe a 
rather shallow minimum in 9 starting around 9 = 1 and ranging up to 9 = 2 where 
the relative errors of the renormalization constants for both operators Ou and O44 
are about the same and comparable to AZ/Z of Zq 12 (9 = 2ir). The case of Zq A4 
is hence an example where using 9 instead of standard lattice momenta can help 
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substantially. From Fig. 0] the value of 9 to be chosen for a practical simulation can 
not be clearly identified. However, taken our analysis of the convergence behavior 
of perturbation theory, the cutoff effects as computed in 1-loop perturbation theory 
and the investigation of the signal to noise ratio together, we conclude that 9 = 1 
should be a reasonable choice. 



3 Non-perturbative step scaling functions 

Choosing the lattice regularization discussed above, is, of course, motivated by the 
advantage to perform numerical simulations. The aim of such simulations is to com- 
pute the lattice step scaling functions of eq. non-perturbatively at non- vanishing 
values of the lattice spacing and to perform the -well-defined- continuum limit of 
eq. (J33J). Besides the lattice step scaling function we will also consider 
constructed in a completely analogous way from the definition of the renormaliza- 
tion factors in eq. (fT5j) . In addition, we will compute the step scaling function of fi, 
eq- m, 

ivm] 

M c V / ACM) 

where c = J f[°\sL, 9)/ J f[°\L, 9) is the corresponding ratio at tree-level. The 
lattice version £ f 1 (a/L; 9) of this step scaling function is again defined in the obvious 
way. To complete the definition of the step scaling functions, certain choices of the 
values for 9, s, x^/L and the quark mass m q have to be made. We fix these in the 
following to 

x 4 /L = 1/2, s = 2, 9 = 1, m q = . (28) 

In order to justify our choice of x 4 /L = 1/2 we show a comparison of the step 
scaling function E^ 44 both at x 4 /L = 1/4 and x 4 /L = 1/2 in fig. 03 Strong lattice 
artifacts are visible for x 4 /L = 1/4. Approaching the continuum limit this effect is 
considerably reduced, but it is clear that E^ 44 is much better behaved as a function 
of a/L when x 4 /L = 1/2 is chosen, which allows for a reliable (linear) continuum 
extrapolation. With the parameter choices in eq. (|28jh we will in the following 
only give the arguments of Zq and oz Q when they differ from these values. In 
particular, with the choice of s = 2 the step scaling function becomes a function 
of only one scale, i.e. <Jz {ni/ ^2) — * ^z {^) an d correspondingly the lattice step 
scaling function will only depend on the scale \x (given by a fixed value of g(L)) and 
a/L. 

As already discussed above, the continuum limit of S^ (^ 2 (L), a/L) is to be 
taken at a fixed value of the renormalization scale \x = 1/L in order to obtain the 
desired continuum step scaling function cr Zo (g 2 (L)). Fixing the scale is realized by 
fixing the renormalized coupling g 2 {L), which in turn can be achieved by changing 
the value of the bare coupling as the value of the lattice spacing (and correspondingly 
a/L) is decreased. Using the SF scheme, the matching points of the bare coupling g 
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Figure 5: Lattice step scaling functions £z 44 (a/L) at £4 = L/A and £4 = L/2 at 
g 2 = 1.8811. The value of 9 is taken from our choices in eq. 



to keep g 2 {L) fixed when the lattice spacing is sent to zero can be found in ^H]- Of 
course, the values of T*z {g 2 {L), a/L) themselves have to be computed in a numerical 
simulation on lattices with increasing number of points N 4 = (L/a) 4 as a — > 0. 

As in our previous work, we will choose two different lattice discretizations, 
standard Wilson fermions and non-perturbatively 0(a)-improved (clover) fermions 
to compute the step scaling function T, Zo (g 2 (L),a/L). Since we want to employ a 
massless renormalization scheme and have to stay therefore at zero quark, we take 
from [IB] the values of k c for clover fermions to stay in the massless limit. For the 
case of Wilson fermions, we determined k c ourselves (see tabic EJ. 

Simulations have been performed on APEmille machines using even-sized lattices 
ranging from 6 4 to 32 4 . We recall that for both actions we used the unimproved 
operator: as a consequence, the results at finite values of the lattice spacing are 
affected by lattice artefacts of O(a), whose precise form depends upon the choice of 
the lattice action. 

For the inversion of the Dirac operator we used the implementation ^J] of the 
SSOR-preconditioned BiCGStab inverter [20J. The gauge fields were generated with 
a hybrid of Cabibbo-Marinari heat-bath and over-relaxation updates: we decided 
to keep the number of over-relaxation steps between a single heat-bath one pro- 
portional to L/a (in practice N or = |L/a). A full update sweep is defined as one 
heat-bath step followed by N OT steps: in order to statistically decorrelate the gauge 
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9sf 


L/a 


P 


K c [Wilson] 


K c [Clover] 


0.8873 


6 


10.7503 


0.134696(7) 


0.130591(4) 




8 


11.0000 


0.134548(6) 


0.130439(3) 




12 


11.3384 


0.134277(5) 


0.130251(2) 




16 


11.5736 


0.134068(6) 


0.130125(2) 


1.0989 


6 


9.5030 


0.136520(5) 


0.131514(5) 




8 


9.7500 


0.136310(3) 


0.131312(3) 




12 


10.0577 


0.135949(4) 


0.131079(3) 




16 


10.3419 


0.135572(4) 


0.130876(2) 


1.3293 


6 


8.6129 


0.138346(6) 


0.132380(6) 




8 


8.8500 


0.138057(4) 


0.132140(5) 




12 


9.1859 


0.137503(2) 


0.131814(3) 




16 


9.4381 


0.137061(4) 


0.131589(2) 


1.5533 


6 


7.9993 


0.140003(11) 


0.133118(7) 




8 


8.2500 


0.139588(8) 


0.132821(5) 




12 


8.5985 


0.138847(6) 


0.132427(3) 




16 


8.8323 


0.138339(7) 


0.132169(3) 


1.8811 


6 


7.4082 


0.142145(11) 


0.133961(8) 




8 


7.6547 


0.141572(9) 


0.133632(6) 




12 


7.9993 


0.140597(6) 


0.133159(4) 




16 


8.2415 


0.139900(6) 


0.132847(3) 


2.1000 


6 


7.1214 


0.143416(11) 


0.134423(9) 




8 


7.3632 


0.142749(9) 


0.134088(6) 




12 


7.6985 


0.141657(6) 


0.133599(4) 




16 


7.9560 


0.140817(7) 


0.133229(3) 


2.4484 


6 


6.7807 


0.145286(11) 


0.134994(11) 




8 


7.0197 


0.144454(7) 


0.134639(7) 




12 


7.3551 


0.143113(6) 


0.134141(5) 




16 


7.6101 


0.142107(6) 


0.133729(4) 


2.7700 


6 


6.5512 


0.146825(11) 


0.135327(12) 




8 


6.7860 


0.145859(7) 


0.135056(8) 




12 


7.1190 


0.144299(8) 


0.134513(5) 




16 


7.3686 


0.143113(7) 


0.134114(3) 


3.4800 


6 


6.2204 


0.149685(15) 


0.135470(15) 




8 


6.4527 


0.148391(9) 


0.135543(9) 




12 


6.7750 


0.146408(7) 


0.135121(5) 




16 


7.0203 


0.145025(8) 


0.134707(4) 



Table 2: Simulation Parameters 
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Figure 6: Lattice step scaling functions S/ 1 (^ 2 (L), a/L) and Ez 44 (g 2 (L), a/L) at 
g 2 = 1.8811. The values of # and X4/L are taken from our choices in eq. 



configurations used to take measurements of our observables, we allowed for a num- 
ber of full sweeps between measurements ranging from 20 (lattices up to L/a = 16) 
to 70 {L/a = 32). 

As usual in this kind of finite-size-step simulations, the signal over noise ratio 
deteriorates for larger values of g 2 , in the fully non-perturbative regime. In order 
to maintain constant, up to a certain degree, the relative statistical errors on our 
observables, we are obliged to increase the statistics. For example, in the clover case 
and on L/a = 32 lattices, the number of measurements has to grow from around 300 
at the smallest values of g 2 up to order 700 at the largest value of the renormalized 
coupling. Accordingly, for lattices with a smaller number of space-time points, the 
number of measurements grows from 500 to 2000. The pattern for the Wilson case 
is similar, albeit with an overall smaller statistical sample. 

In Fig. El we show separately the lattice step scaling functions E^ 44 (g 2 (L), a/L) 
and Y,f 1 (g 2 (L),a/L) at four values of a/L for g 2 fixed. Both step scaling functions 
show rather large lattice artefacts. For E f 1 this behavior is in sharp contrast to the 
case of choosing 9 = 0, where the a/L dependence of E^ is basically flat, at least 
in the 0(a)-improved theory [31 E]. This flat behavior led us, in our earlier work, 
to the conclusion to perform the continuum limit of S^ 12 and E/ x separately and 
to compute the desired continuum step scaling function a Zl2 = a z 12 /°'h a fter the 
continuum limit of E and E f 1 had been performed. 
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Figure 7: Examples of the continuum limit of step scaling functions for both Wilson 
and clover fermions. The (constrained) fit performed is linear, taking the 3 data 
points with smallest values of a/L. We also indicate the values of the running 
coupling that correspond to these step scaling functions. 



The example of the step scaling functions, shown in Fig. indicates that for 9 
non- vanishing the strategy might be different. Indeed, we found in the analysis of 
our data -for 9 = 1- that in the ratios 

Z 2 Jg 2 (L), a/L)/E fl {f{L), a/L) (29) 

and 

Z 2l2 (g 2 (L),a/L)/Z fl (f(L),a/L) (30) 

the lattice artefacts essentially cancel. This cancellation happens for clover as well 
as for Wilson fermions. We decided therefore that for our choice of 9 = 1 to first 
compute 

E z {f(L),a/L) = E 2 {f{L),a/L)/E fl (f(L),a/L) (31) 

at a given value of a/L and then perform the continuum limit. 

In Fig. [7| we show examples of the continuum extrapolations of £e> 44 for both 
discretization, Wilson and clover fermions. It can be observed that in both cases 
the continuum extrapolation is linear in the lattice spacing when the data points 
with the three smallest values of a/L are taken. In the plot we show already a 
constraint fit, demanding that both set of lattice results extrapolate to the same 
value of the continuum step scaling function. We obtain a very similar plot for 
Ee) 12 . The numerical results of the continuum extrapolation as well as for the four 
lattice spacings are summarized in table 01 for the On and in table |U for the O44 
operator. 
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Figure 8: Continuum step scaling function cr^ 44 . The solid line represents the 
parameterization eq. (|32jh the dotted line represents the the ldoop, the dashed line 
the 24oop predictions from perturbation theory. 

One subtlety, we want to mention, is that even in the continuum limit at fixed 
scale, the step scaling functions of O44 and O12 do not assume the same value. The 
reason is that both operators belong to two different representations and in the here 
adopted finite volume renormalization scheme they are therefore also part of the 
precise definition of the scheme. It is hence not possible, to constrain the fits even 
more by demanding that £e> 44 and £e> 12 converge to the same value of the continuum 
step scaling function. 

4 Running and invariant step scaling function 

At this stage, we can leave the lattice and discuss continuum quantities only. The 
only reminder that we have performed a lattice calculation is that we will stay in 
the somewhat unusual, finite volume SF renormalization scheme. The simulations 
described in the previous section provide us with the continuum step scaling func- 
tions cr^ 44 and o~z 12 at a number of values for the running coupling 0.8 < g 2 < 3.5, 
corresponding to a wide range of scales, 500 MeV < \i < 100 GeV. 

In Fig. |H1 we show the continuum step scaling function o~z i4 as a function of the 
running coupling (and hence the scale). The dotted curve is the ldoop perturbative 
analysis for the step scaling function, whereas the dashed curve is a 2- loop result, 
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Figure 9: Continuum step scaling function <Jz 12 - The solid line represents the 
parameterization eq. (|32|). the dotted line represents the the 1-loop, the dashed line 
the 2-loop predictions from perturbation theory. We also show the values of az l2 for 
6 = 2n (open symbols) corresponding to a minimal physical momentum. 

taking the 2-loop anomalous dimension in the SF scheme from its relation to the 
2-loop anomalous dimension in the MS-scheme discussed in section 2.1. The solid 
line is a fit to the data according to the formula 



In eq. (J3*2*)l we have taken the known 1-loop (70) and fitted only the coefficients ci, c 2 
which we find to be C\ = —0.0334(50), c 2 = 0.0041(18). It can be seen from the 
figure that 1-loop perturbation theory is not a good description of the data at all. 
Even 2-loop perturbation theory represents the data only up to g 2 = 2, or so. 

Fig.EJis the same as fig.|Hlfor the operator 12 . In this case we find the coefficients 
ci = -0.0343(41), c 2 = 0.0049(14) for the parameterization of eq. ([52]). In addition, 
we show the values of the step scaling function a Zi2 as obtained in our earlier work, 
i.e. using 8 = and a momentum p = 27r /L (open symbols). Clearly, the scale 
dependence of cr Zl2 (9 = 1) is much stronger than for & Zl2 - 

The final result of this paper is the (ultra-violet) invariant step scaling function 
for the operators C44 and Ou as they are needed to describe the running of the 
physical matrix element in other schemes than the SF one. We adopt here the 



0-Z44 = l-7o ln(2)#o + ci#o + c 2 #o • 



(32) 
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definition of the invariant step scaling function as given in [Sj: 



©Sv^M = 9 2 (L )) ■ (f(L)r^ 2b0 exp j - dg M - J± 

' (33) 

The (3- and 7-functions in eq. f|33|) will be taken to a given order in perturbation 
theory in the SF scheme. A scheme dependence is still manifest in ©inv through 
the appearance of the infra-red scale (/x = -^o 1 )- We thus add a superscript to 
indicate that the invariant step scaling function is computed within the Schrodinger 
functional scheme. 

For large enough scales in the ultra-violet, it is expected that the running of 
the step scaling functions, as computed non-perturbatively above, will match the 
perturbative running and hence that ©inv w iU become independent from the scale. 
In order to test at which scales ©inv becomes constant, we have to evolve the step 
scaling functions starting from our most non-perturbative scale fiQ = (l/2L max ), 
implicitly given by g 2 (2L max ) = 3.48. One finds [2T] 2L max = (1.476r ), where 
r ~ 0.5fm. The step scaling functions <j Zo are only computed at certain values of 
g 2 (L). Subsequent values of g 2 (L) do, however, not correspond to a scale change by 
a factor of two as needed. In order to find the precise evolution by steps of two an 
interpolation in g 2 (L) has to be performed using the parameterization of eq. (|32|). 
The corresponding values of the step scaling function can be found in table EJ 

The error of the so obtained values of oz Q is evaluated by a standard error 
propagation taking the correlation of the fit parameters into account through the 
covariance matrix (cov), i.e. 



(A. Zo ) = ^ W + ^ W + ^ • ^ • cov, + . 

(34) 

As an aside we mention that we also have seen that the errors coming from the 
uncertainty in g 2 (L) are not negligible, and they are included in the values given in 
table El In evaluating ©inv i n ec L- we have taken the 3- loop /9-function and the 
2-loop 7-function. 

In Fig. [TUlwe show ©inv as a function of fi/A^ with w 120 MeV b eing the 
A-parameter in the quenched approximation in the SF scheme. For the operator O44 
(full symbols) the invariant step scaling function is constant for, say, /x/Agp > 50 
indicating that contact with perturbation theory can safely be made. In the same 
figure we also show the value of the invariant step scaling function for Ou which 
shows a very similar behavior. The fact that both invariant step scaling functions 
assume different values indicates again that the two different operators define two 
different renormalization schemes, at least in a finite volume scheme like the SF. 
The difference will only disappear in the physical renormalization group invariant 
matrix element. 
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Figure 10: The RG invariant step scaling function for the operator O44 (full symbols) 
and the operator C 12 (open symbols). 



In order to obtain finally the values of the invariant step scaling functions, we 
have taken the values at the corresponding largest scale: 

©invXM = 0-242(8), ©g&f£> ) = 0.221(9) . (35) 

As a comparison we quote the value of ©rNV S a L2 (/ i °) = 1-1 1(4) as obtained at 9 = 
and p = 2tt/L in our previous work. 



5 Conclusion 

In this paper we have demonstrated how the generalized boundary conditions ad- 
vocated in ref. pQ can to utilized to define continuous external momenta as often 
needed in the lattice simulations when operators with derivatives are considered. We 
investigated the particular example of a scale dependent renormalization constant 
for renormalizing twist-2 non-singlet quark operators corresponding to moments of 
parton distribution functions. 
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Although we demonstrated the usefulness for such momentum definitions for a 
renormalization constant, or more specific its step scaling function, it is clear that 
in principle the same method also applies to matrix elements themselves, as long as 
no real physical momentum transfer is involved. The feasibility of this approach has 
to be tested, however, in a real numerical benchmark simulation. 

For the particular example at hand, we could determine the renormalization con- 
stant for a certain lattice representation of a continuum operator that determines 
the average (quark) momentum in a hadron. Such a calculation would have been 
very difficult by using standard, quantized lattice momenta as usually taken. As 
a result, we could determine the RG invariant step scaling function for two differ- 
ent lattice representations of the same continuum operator. In this way, we could 
eliminate a small systematic uncertainty that plagued our earlier work. 

The idea of the present work was to show the practicability and usefulness of 
the generalized momentum definition on a practical example. The values of the 
RG invariant step scaling we have computed are essential ingredients for the RG 
invariant matrix elements which allow a comparison with experimental data or global 
fits. Results for these physical matrix elements, in which we are finally interested 
in, will be reported in a forthcoming publication. 
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9sf 


a/L 


Ez 12 [Wi/son] 


£z 12 [Clover] 


a z 12 


X 2 /dof 




1/6 


0.9937(25) 


0.9682 


;60) 






0.8873 


1/8 


0.9787(34) 


0.9616 


;59) 


0.9475(96) 


0.32 




1/12 


0.9640(55) 


0.9571 


;58) 








1/16 


0.9689(94) 


0.9555 


71) 








1/6 


0.9899(39) 


0.9540 


;58) 






1.0989 


1/8 


0.9727(55) 


0.9591 


;62) 


0.9211(105) 


0.23 




1/12 


0.9608(67) 


0.9446 


;6i) 








1/16 


0.9480(91) 


0.9368 


74) 








1/6 


0.9771(42) 


0.9360 


'69) 






1.3293 


1/8 


0.9680(61) 


0.9365 


;67) 


0.8905(120) 


0.79 




1/12 


0.9339(81) 


0.9305 


;66) 








1/16 


0.9230(117) 


0.9157 


;8i) 








1/6 


0.9716(48) 


0.9153 


;64) 






1.5553 


1/8 


0.9545(77) 


0.9102 


[67) 


0.8896(129) 


0.26 




1/12 


0.9287(92) 


0.9068 


[63) 








1/16 


0.9133(132) 


0.9031 


[86) 








1/6 


0.9607(68) 


0.8984 


75) 






1.8811 


1/8 


0.9259(104) 


0.8778 


70) 


0.8522(144) 


1.07 




1/12 


0.9060(103) 


0.8781 


72) 








1/16 


0.8682(153) 


0.8640 


;93) 








1/6 


0.9627(68) 


0.8700 


71) 






2.1000 


1/8 


0.9107(74) 


0.8592 


70) 


0.8419(131) 


1.01 




1/12 


0.8981(101) 


0.8558 


72) 








1/16 


0.8827(102) 


0.8368 


;99) 








1/6 


0.9467(67) 


0.8217 


75) 






2.4484 


1/8 


0.8811(85) 


0.8223 


77) 


0.8064(133) 


0.78 




1/12 


0.8589(136) 


0.8053 


77) 








1/16 


0.8519(85) 


0.8116 


107) 








1/6 


0.9309(75) 


0.7864 


75) 






2.7700 


1/8 


0.8759(94) 


0.7838 


78) 


0.7395(157) 


0.77 




1/12 


0.8385(102) 


0.7593 


75) 








1/16 


0.8225(152) 


0.7583 


118) 








1/6 


0.9063(50) 


0.7176 


'83) 






3.4800 


1/8 


0.8203(110) 


0.6899 


;86) 


0.6793(187) 


0.56 




1/12 


0.7846(144) 


0.6761 


;89) 








1/16 


0.7604(191) 


0.6851 


136) 







Table 3: Results for the lattice step scaling function Ti Zl2 and the combined contin- 
uum extrapolation a Zl2 of Wilson and Clover data with the three smallest lattice 
spacings. 
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9sf 


a/L 


ZzjWilson] 


S Z44 [Clover] 


a z 44 


X 2 /dof 




1/6 


0.9776(30) 


0.9596(68) 






0.8873 


1/8 


0.9677(44) 


0.9513(68) 


0.9489(115) 


0.49 




1/12 


0.9671(65) 


0.9514(68) 








1/16 


0.9601(125) 


0.9425(81) 








1/6 


0.9577(49) 


0.9375(66) 






1.0989 


1/8 


0.9560(69) 


0.9421(71) 


0.9234(128) 


0.50 




1/12 


0.9505(85) 


0.9420(71) 








1/16 


0.9330(123) 


0.9290(86) 








1/6 


0.9402(54) 


0.9184(77) 






1.3293 


1/8 


0.9416(81) 


0.9275(77) 


0.8943(145) 


1.40 




1/12 


0.9250(104) 


0.9322(78) 








1/16 


0.9118(153) 


0.9037(93) 








1/6 


0.9062(63) 


0.8908(71) 






1.5553 


1/8 


0.9163(100) 


0.8977(75) 


0.8900(154) 


1.60 




1/12 


0.9140(117) 


0.9084(74) 








1/16 


0.8834(163) 


0.8873(103) 








1/6 


0.8813(88) 


0.8777(84) 






1.8811 


1/8 


0.8723(140) 


0.8569(80) 


0.8407(178) 


0.54 




1/12 


0.8723(139) 


0.8577(84) 








1/16 


0.8434(199) 


0.8439(117) 








1/6 


0.8698(88) 


0.8354(80) 






2.1000 


1/8 


0.8467(97) 


0.8399(82) 


0.8373(157) 


1.75 




1/12 


0.8545(133) 


0.8425(84) 








1/16 


0.8590(129) 


0.8181(112) 








1/6 


0.8197(88) 


0.8016(89) 






2.4484 


1/8 


0.7935(119) 


0.7885(91) 


0.7944(169) 


0.55 




1/12 


0.7942(186) 


0.7921(94) 








1/16 


0.7823(115) 


0.8036(127) 








1/6 


0.7589(103) 


0.7494(88) 






2.7700 


1/8 


0.7606(129) 


0.7556(90) 


0.7376(197) 


1.22 




1/12 


0.7736(143) 


0.7347(92) 








1/16 


0.7542(203) 


0.7512(143) 








1/6 


0.6739(66) 


0.6528(99) 






3.4800 


1/8 


0.6376(153) 


0.6469(104) 


0.6525(241) 


0.34 




1/12 


0.6500(193) 


0.6385(110) 








1/16 


0.6522(256) 


0.6563(178) 







Table 4: Results for the lattice step scaling function E^ 44 and the combined contin- 
uum extrapolation a Zii of Wilson and Clover data with the three smallest lattice 
spacings. 
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i>SF 


(J7,„ 
^12 


<7 7 


3.480 


0.6826(169) 


0.6579(216) 


2.454(18) 


0.7893(83) 


0.7822(100) 


1.918(18) 


0.8486(76) 


0.8458(90) 


1.584(18) 


0.8840(67) 


0.8828(79) 


1.353(18) 


0.9071(58) 


0.9066(68) 


1.184(17) 


0.9231(50) 


0.9229(58) 


1.053(15) 


0.9348(42) 


0.9348(49) 


0.950(14) 


0.9436(37) 


0.9437(42) 


0.865(13) 


0.9505(32) 


0.9506(37) 



Table 5: Interpolated values of the step scaling function a Zl2 and a^ 44 for the scale 
change by a factor two. 



Appendix A 

In this appendix we recall some basic properties of the renormalization group func- 
tions. For small couplings the /5-function and the anomalous dimension (7-function), 
defined by 

d 

J (36) 
have asymptotic expansions of the form 

00 

k Z° (37) 

1(9) g ~°-g 2 J2^ 2k - 

k=0 

One finds that b , b\ and 7 are the same in all the schemes (these are the "univer- 
sal" coefficients), while all other coefficients are scheme dependent. The universal 
coefficients are given by (with C F — (AT? — 1)/2N C and N c the number of colors) 

bo = (yiV c -jjiV f }(47r)- 2 , (38) 
1 6 

7o = -jC F (4n)- 2 , (39) 
34 1 3 

h = {-iV c 2 -(-iV c -iV- 1 )iV f }(47r)- 4 . (40) 
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Any two mass independent renormalization schemes can be related by a scale change 
and a finite parameter renormalization of the form 

// = Cfi, c> 0, (41) 

g' = gyfx^l (42) 

Z' = Z[AZ (g)}, (43) 

where c is just a change of scale between the 2 schemes and one could obviously 
choose also c = 1. X„ and AZo are expanded according to 



xM + ( 44 ) 



AZ (^)^°l + ^A4 fc) ^. (45) 
fc=i 

The invariance of a physical observable under such a change of parameters, gives a 
relation between the renormalization group functions, (3 and 7, in the 2 schemes. In 
particular we have 

7l = 7 p + 2&0 AZ ( g - 70 * g (1) • (46) 

From ref. [22] we have 

* g (1) = -7^(ci,o + ci,iiVf), (47) 

6 47T 

with 

ci, = 1.25563(4), ci,i = 0.039863(2). (48) 
From the perturbative results in [TBj (section 2.1) it is possible to obtain AZq\ 
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